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^^ , Abstract. Masses ol black holes in nearby galactic nuclei can be measured in a 

variety of ways, using stellar and gaseous kinematics. Reliable black hole masses 
are known for several dozen objects, so that demographic questions can start to be 
f^ l- addressed with some confidence. Prospects for the near future are discussed briefly. 
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1 Introduction 

Active galaxies and quasars are powered by physical processes in an accretion 



disk surrounding a massive black hole p7 54 1 . The observed number of active 
galaxies increases towards high redshift z, so that many 'normal' galaxies 
^k.' I must have been active in the past p9| |. This implies that inactive massive 

^S| ■ central black holes must lurk in the nuclei of nearby normal galaxies. In the 

0^ I past decade, much work has been done to measure the masses of these black 

^p . holes, to establish the relation between black hole mass and the global/nuclear 

^^ ' properties of the host galaxy, and to understand the role these objects play 

in driving internal dynamical evolution |{47| , |64| . 

A black hole of mass A/bh in a galactic nucleus dominates the gravi- 
tational potential inside the so-called radius of influence which is usually 
O ' defined as tbh — GA/bh/c^, where G is the gravitational constant, and a is 



the characteristic velocity dispersion in the host galaxy. In physical units 

„,/ Mbh \/100km/sx2 

For a galaxy at distance D, tbh corresponds to an angular size 

/ A/bh \ /100km/s\2/lMpc 



Inside the sphere of influence, the black hole generates a central cusp in the 
density distribution of the stars. The resulting density and luminosity oc r~'^ 
with 3/2 < 7 < 9/4 [^. The typical velocities scale ex r~^/^ with radius. 

While many early-type galaxies have cusps in their central luminosity 
profiles I9p8,p3[ , these by themselves are not proof of the presence of a black 



hole, as other processes can generate such cusps |34 . For this reason most 
studies concentrate on obtaining spectroscopy at the smallest angular scales, 
to find evidence for the expected Keplerian rise in the velocities of stars 
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Fig. 1. Tlie radius of influence ran = GMbh/o"^ of central blaclc holes with well- 
determined masses versus the distance of the host galaxy. The masses were derived 
from stellar kinematics (stars), gas kinematics (dots) and VLBI measurements of 
masers (triangles). The solid line corresponds to an angular size Obh ~ I'.'O, and 
the dotted line corresponds to O'.'l. The error bars represent the quoted uncertainty 
in AfBH @|jI|,|I|,0,§j2i[g,§,§J38m|,||,||. 



and gas. This requires high spatial resolution. For example, the 3 x IO^Mq 
black hole in the nucleus of the galaxy M87 [^ in the Virgo cluster has 
^BH ~ 1''5. The general approach is to determine the luminous mass in stars 
from the observed surface brightness distribution, and to compare this with 
the dynamical mass derived from kinematic measurements of stars or gas 
[ ^2|j55[ | . If one can show that the mass density inside a certain radius is larger 
than anything that can be produced by normal dynamical processes, then 
the object is considered to be a black hole. In a few cases it is possible to 
find direct evidence for the presence of a relativistic object (§|[). 

Figure |l| shows tbh as defined in eq. (1^) for the best published black hole 
mass determinations versus the distance to the host galaxy. Lines of constant 
angular resolution ^ros run diagonal. The early determinations clustered near 

O'.'l. 



I'.'O l56 



but HST has pushed this to Ores ~ O'.'l. Detailed modeling 
has shown that, depending on the internal dynamical structure of the host 
galaxy, the effects of the black hole often are visible only inside projected radii 
that are significantly smaller than tbh (e.g., eq. (4.2) in jS^). This suggests 
that measured masses corresponding to 0bh ^ ^ics should be treated with 
caution, as they are likely to be overestimates (see i^). At present, only VLBI 
measurements can probe the regime 0res < O'.'l. 
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2 Stellar dynamical modeling 

The nearest galactic nucleus is the center of our own Galaxy. Despite the large 
foreground extinction, it is possible to resolve individual stars in the Galactic 
center in the infrared, and to measure not only their radial velocities, but also 
their proper motions |15| , p8| , and accelerations [ P9| ! Dynamical modeling of 
this remarkable data provides unequivocal proof that our Galaxy contains a 
central black hole of nearly three million solar masses. 

The dynamics of the nuclei of nearby early-type galaxies can be probed 
with stellar absorption- line spectroscopy of the integrated light. This gener- 
ally requires long exposure times. The orbital structure in these systems is 
rich |^,g4|, so a true inward increase of the mass-to- light ratio M/L must 
be distinguished from radial variation of the velocity anisotropy. This can be 
done by measuring the shape of the line-of-sight velocity distribution ||2^,^ . 
The orbital structure is related to the intrinsic shape of the galaxy. This can 
be constrained by measurements along multiple position angles P,p4|, or, 
even better, by integral-field spectroscopy (§||). 

Determination of the black hole mass and the orbital structure requires 
construction of dynamical models. There has been a steady increase in the 
sophistication of model construction in the past decade. Early isotropic spher- 
ical models were replaced by anisotropic spheres, and then by axisymmctric 
models with a special orbital structure (phase-space distribution function 
/ = f{E^ Lz) where E is the orbital energy and L^ is the angular monrcntum 
component parallel to the symmetry axis, e.g., [p2[). More recently axisym- 
mctric models with the full range of possible anisotropics, and multiple com- 
ponents have been used. The first such study was done for M32, and included 
ground-based data along four position angles and eight FDS pointings |Q. 
The model was constructed by a version of Schwarzschild's ||57| numerical 
orbit-superposition method, which fits the surface brightness distribution as 
well as all kinematic observables [|l3[ . 

Another example of this approach is provided by the E7/S0 galaxy NGC 
4342 in Virgo. This is a low-luminosity object, seen nearly edge-on, with a 
prominent nuclear stellar disk. Cretton & van den Bosch |Q used ground- 
based major-axis kinematics from the WHT and five FOS pointings 0, and 
compared these with general axisymmctric models containing a spheroid, a 
stellar disk, a nuclear disk and a black hole. Free parameters in the modeling 
were M-qu and the stellar mass-to-light ratio T. Contours of constant x^ in 
the (T, MBH)-plane show that the best fit is obtained for A/bh -^ 3 x lO^M©, 
but with a significant uncertainty (Figure ||). 

To date, black hole masses have been derived in this way for about fif- 
teen objects. The full analysis has been published only for M32 Q, N4342 
[^, and NGC3379 [|4), all based on ground-based and FOS data. Black hole 
masses based on three-integral axisymmctric modeling of STIS absorption- 
line spectroscopy have been reported for a dozen objects by Gebhardt's group 
[p5|, but the data and the models have not yet been published. 
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Fig. 2. Tiie black iiole in NGC 4342. Contours indicate the goodness-of-fit x^ to 
the observed photometry and kinematics as a function of black hole mass Mbh and 
T/, the 7-band stellar mass-to-light ratio. Dots indicate dynamical models that were 
constructed. The asterisk denotes the best fit. The first three contours surrounding 
it are the 68.3%, 95.4% and 99.7% confidence levels, while the subsequent contours 
correspond to a factor of two increase in x^ P^l • 



While in some cases axisymmctric models are consistent with the data, in 
others they clearly are inadequate (§^. The nucleus of M31 has been known 
to be asymmetric since 1974 [3^. Observations with the integral-field spec- 
trograph TIGER revealed that the black hole resides in the secondary peak in 
the brightness distribution p], and led to models with an asymmetric distri- 
bution of stars |£^. Measurements with OASIS on the CFHT [|| have shown 
that HST spectroscopy with FOC ^^, FOS and STIS along the apparent sym- 
metry axis of the eccentric structure has in fact missed the kinematic major 
axis. A comprehensive non- axisymmctric model that fits this data will teach 
us much about the structure and formation of this nearby nucleus. 

Some nearby galactic nuclei are shrouded in dust, and their internal kine- 
matics is best probed at longer wavelengths. This is now possible through 
the availability of near-IR spectrographs, which employ the CO bandhead 
at 2.3/x to derive the stellar kinematics. Anders Q] used the MPE-built 3D 
integral-field unit to show that the derived nuclear kinematics in the largely 
dust-free nearby galaxy NGC 3115 agrees with the kinematics measured at 
shorter wavelengths W& , demonstrating that this approach works. 3D obser- 
vations of the luminous merger remnant NGC 1316 (Fornax A) confirmed the 
central a of w230 km/s obtained at short wavelengths from the ground and 
with FOS [|||, suggesting Mbh ^ lO^Mg Q. Instruments such as SINFONI 
will probe nearby dusty nuclei with a resolution similar to HST Pq|. 
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3 Gas kinematics 



Optical emission lines. The nuclei of active early-type galaxies, and those of 
most spirals, contain extended optical emission-line gas. Its kinematics can 
be used to constrain the central mass distribution. In this case the exposure 
times can be relatively short, but care must be taken to model possible non- 
circular motions and the effects of turbulence in the gas disks. 

High spatial resolution emission- line kinematics (FOC, FOS and STIS), to- 
gether with careful modeling, has been published for six cases: M87, NGC 
4261, NGC 4374, NGC 6251, NGC 7052, IC 1459 |J20l|23,||,||J6| . These 
very luminous early-type galaxies cover a modest range in total luminosity, 
but the derived black hole masses vary by a factor of ten. 

An example of this approach is provided by the E3 galaxy IC 1459. This 
giant elliptical hosts a compact nuclear radio source, has a counter-rotating 
stellar core ('^ 10"), a shallow cusp in the luminosity profile, and a blue cen- 
tral point source |^,^. FOS kinematics of the emission-line gas associated 
with the nuclear dust-lane reveals a disk in rapid rotation (Figure ra), with 
significant velocity dispersion. Detailed models for the gas motions which in- 
clude the effects of turbulence show that AIbh ~ 2.5 x IO^Mq |p3|, which 
is a factor of ten smaller than the mass suggested by axisymmetric f{E, L^) 
models of the ground-based stellar kinematics along the major axis. A general 
stellar dynamical model that incorporates major-axis STIS spectroscopy and 
ground-based spectroscopy along four position angles is under construction. 
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Fig. 3. Emission-line gas kinematics of the nucleus of IC1459, derived from three 
FOS pointings with the C/.'l (left) and three with the 0'.'25 apertures (right). Rota- 
tion velocities V were derived from Hq-|-[NII] (open circles) and H/3 measurements 
(dots). The heavy solid line is the prediction of a model with Mbh = 1.0 x 10 Mq. 
Dotted and dashed curves are for Mbh = and 7 x 10* Mq, respectively ||63|. 
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Kinematics of emission-line gas will provide many black hole masses in the 
near future. Ongoing HST programs include Ha emission-line spectroscopy 
with STIS for 21 radio-loud ellipticals (PI Baum) and 54 Sb/Sc spirals (PI 
Axon). These studies generally use three parallel slits to measure possible 
deviations from simple circular motion. 

Masers. Nuclear maser emission can be measured with VLBI techniques in a 
few spiral galaxies. This achieves the highest spatial resolution to date, but 
is possible only when the circumnuclear disk is nearly edge-on. A search of 
«700 nuclei revealed maser emission in 22 cases, and disk-like kinematics in 
six of these fpOf- The best black hole mass determinations are for NGC 1068, 



NGC 4258, and NGC 4945 [BO 49 



X rays. Recently, it has become possible to measure the profile of the Fe Ka 
line at 6.4 keV in the nuclei of nearby Seyfert galaxies |pl|. The width of 
this line approaches wlO^ km/s, which is direct evidence that the emitting 
gas must be near the Schwarzschild radius of a relativistic object. The mass 
of the black hole, and possibly its spin, can be derived from the detailed 
shape of the line profile [p9p^. The latest generation of X-ray telescopes, in 
particular CHANDRA and XMM, will extend this work to more objects. 

Reverberation mapping. The observed time-variation of broad H/3 emission 
lines from Seyfert nuclei can be used to derive the radius of the broad-line 
region by means of the light travel-time argument [^ . In combination with 
simple kinematic models for the motion of the broad- line clouds, this provides 
an estimate of the mass of the central object responsible for these motions. 
Early mass determinations appeared to be systematically lower than those 
obtained by other means |33] , but this discrepancy has now disappeared ||26| . 



4 Demographics 

For the past five years, our understanding of black hole demography was 
summarized in a diagram of black hole mass Mbh versus absolute bulge lu- 
minosity Lbuigc, where the 'bulge' was taken to be the entire galaxy for an 
elliptical or lenticular, or the actual bulge in case of a spiral [p5U5q]. Early 
work based on simple axisymmetric Jeans models and ground-based kine- 
matics [p9| was interpreted as evidence for a tight correlation of A/bh with 
Lhu\ge, and hence with the total mass of the bulge. Much effort was spent 
in trying to reproduce this correlation in models of galaxy formation, and in 
relating it to the energy production in quasars |55|] . 

More appropriate dynamical models combined with HST kinematics have 
generally resulted in a downward revision of the earlier masses |2aJ33,Ert. 
Figure Ha is the resulting demography diagram for the best determinations. 
It shows a rough correlation, but with a large range in Mbh at fixed Lb- The 
masses are consistent with various black hole formation scenarios pc| , |6l[|, and 
they agree with the quasar light prediction for reasonable efficiencies ]55|,p9{ . 
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Fig. 4. Demography of black holes. Left: A/bh versus total absolute luminosity 
Lb, bulge of the host bulge or spheroid. Coding of points is as in Fig. m. The dashed 
line is the correlation based on f(E, Lz) models [M. The solid line is the prediction 
for adiabatic growth models ^0|). Recently published measurements generally give 
lower masses, indicating mild radial anisotropy. They display a significant scatter 
at fixed Lb, bulge- There is an observational bias against detecting small black hole 
masses in luminous galaxies. Right: Mbh versus velocity dispersion o of the host 
bulge or spheroid. The solid line is the relation proposed in psl. 



However, there is an observational bias against detecting small black holes in 
large galaxies, i.e., for cases where the radius of influence is simply smaller 
than the spatial resolution currently achievable. 

Absolute galaxy luminosity correlates with velocity dispersion |17|| , so 
Mbh is expected to correlate with a, measured outside the radius of influence 
of the black hole. As illustrated in Figure [ip, the scatter in this correlation is 
considerably less than in the Mbh versus Lbuige correlation [p2|J25|| . A good fit 
is provided by Mbh = 1.30(±0.36) x lO8M0(cr/2OOkm s"i)4"^0-36)^ where 
a is the velocity dispersion at Re/8, and Re is the effective radius |Q. This 
suggests yet another link between global and nuclear galaxy properties, but it 
is not understood what causes this link, or whether there are other parameters 
involved. Some of the galaxies in the diagram are not axisymmetric, so the 
black hole mass derived from axisymmetric modeling may need revision. 

Substitution of the Mbh versus a correlation in eq. (||) provides the fol- 
lowing estimate of the formal radius of influence of the black hole: 



7BH 



0' 



/.'a( 



100 km/s 



2.7 



IMpc 
D 



(3) 



As noted in §gj, the actual radius of influence may be smaller. Assuming 
that eq. (^) is valid for all galaxies (which is not proven) allows one to use 
a ground-based velocity dispersion measurement to estimate the minimum 
resolution needed to measure a reliable black hole mass. 
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Fig. 5. Integral-field spectroscopy of tlie E6 galaxy NGC 3377. Top, from left to 
right: reconstructed intensity, stellar Mg b index, mean velocity, and velocity dis- 
persion. Bottom, from left to right: gas intensity (OIII A5007) and velocity, and 
stellar mean velocity and velocity dispersion of the nuclear region. The first six 
panels are based on a two- hour exposure with SAURDN. The field shown is 30" x 39". 
The last two panels are based on a 3.5 hour exposure with OASIS, with 0'.'l6 sam- 
pling and a 4" x 4" field. The stars show a striking rotating disk pattern with the 
spin axis misaligned ~ 10° from the photometric minor axis, indicating the galaxy 
is triaxial. The gas also reveals non-axisymmetric structures and motions |lO[|ll|]. 



5 Next steps 



There is strong evidence for the existence of massive black holes in most 
galactic nuclei, based on different measurement approaches. Our understand- 
ing of black hole demographics is still limited, because the number of reliable 
mass measurements is relatively modest. Significant improvement is expected 
through spectroscopic surveys of galactic nuclei with emission-line gas, and 
through further work on stellar absorption-line kinematics. 

To make further progress, high-resolution observations of the nuclei ob- 
tained with STIS onboard HST, and with adaptive optics from the ground 
(e.g., OASIS on the CFHT, SINFONI on the VLT), need to be complemented 
with the wide-field kinematics of the host galaxy, in order to measure its 
intrinsic shape and orbital structure [ppj]. For this reason, the dynamics 
groups of the universities at Lyon, Leiden, and Durham have built the special- 
purpose integral-field spectrograph SAURDN for the WHT, with high through- 
put, and a field of view of 33" x 41" sampled at 0"94 x 0'.'94 [|,||, and are 
using it to observe a representative sample of 80 nearby early-type galaxies. 

Figure S shows integral-field kinematics for the E6 galaxy NGC 3377 
[10[jll| . The SAURDN maps of the stellar and gaseous motions reveal that this 
galaxy is not axisymmetric. The minor axis rotation persists in the high- 
spatial resolution measurements with OASIS. NGC 3377 has a steep-cusped 
central luminosity profile, and the expectation based on general dynamical 
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arguments is that this galaxy should be nearly axisymnietric, at least in the 
inner regions [ p7|j64| . It will be interesting to probe this galaxy at even higher 
spatial resolution, to establish whether the nucleus is axisymmetric and the 
derived black hole mass correct, or whether the non-axisymmetry persists, 
and our understanding of the nuclear dynamics needs modification. Either 
way, the codes for numerical model construction need to be generalized to 
triaxial geometry, and to asymmetric systems such as the nucleus of M31. 

The ongoing systematic programs will reveal the black hole demographics 
as a function of Hubble type, radio properties, intrinsic shapes, and internal 
dynamics. It will also establish the rate of occurrence of gaseous and stellar 
nuclear disks, and of nuclear star clusters, and the importance of black hole 
driven secular dynamical evolution for shaping galaxies. 

It is a pleasure to acknowledge comments and contributions by Nicolas Cret- 
ton, Roeland van der Marel, Gijs Verdoes Kleijn, and the SAURON team. 
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